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ABSTRACT 

This is a second paper of a series of two. In this paper, we directly correct the problem pointed out in 
the first paper of this series, dark energy does not cluster on the scale of galaxy clusters, but the current 
describing the flowing of dark energies outside the clusters is ignored in almost all the previous papers. 
We set up and solve a first order differential equation which describes the evolution of the clusters in a 
back ground universe containing dark energies. From the solution we extract the key parameters of the 
model and find them depending on the equation of state coefficients of dark energies rather non-trivially. 
We then apply the results in Press-Scheter theory and calculate the number density of galaxy clusters and 
its evolutions, we find the observable quantities are strongly affected by the equation of state coefficients 
of dark energies. 

Subject headings: top-hat spherical collapse model, galaxy clusters formation, dark energy 



1. Introduction 

In a previous paper (Zeng & Gao 2005), we point 
out that, in most of the existing literatures about top- 
hat spherical collapse model of galaxy clusters for- 
mation, (Barrow & Saich 1991; Eke, Cole & Frenk 
1996; Wang & Stein hardt 1998; Lokas & Hoffman 
2001; Percival, Miller & Peacock 2000; Weinber & 
Kamionkowski 2002), the dark energy is assumed not 
to cluster on the scale of galaxy clusters. But the cur- 
rent which describes the flowing of dark energy outside 
the clusters is usually ignored, so the discussions in 
these literatures are not self-consistent. In that paper, 
by assuming that dark energy clusters synchronously 
with ordinary matters so that the dark energy current 
flowing outside the clusters does not exist at all, we 
make our discussions self-consistent. However, by so 
doing, we do not correct the problem in the existing 
literatures directly. We only from the contrary indi- 
cate that the effects of the dark energy current may be 
important. 

The purpose of this paper is to directly consider 
the effects of such a current. Just as we state in that 
paper, when we add such a current in the energy- 
momentum tensor, Einstein equation becomes compli- 
cated and the metric of the cluster-inside space-time 
cannot be factorized as usual. So not only Friedman 



equation, but also Raycharduri equation does not fol- 
low as thought by the author of (Wang & Stein hardt 
1998) and (Weinber & Kamionkowski 2002). 

Using energy conservation law, we set up a first or- 
der differential equation to describe the evolution of 
the radius of a cluster and solved it in the subsection 
2.1. We then give numerical as well as formal analyt- 
ical solutions for this equation. In subsection 2.2, we 
extract the key parameters of this model which will 
be used in the application of Press-Scheter theory. In 
subsection 3.1 we derived theoretical formulaes to cal- 
culate the number-density v.s. temperature of galaxy 
clusters using Press-Scheter theory. In subsection 3.2, 
we provide numerical results for the number density 
of galaxy clusters and its evolutions and study the ef- 
fects of dark energy equation of coefficients on this two 
quantities. We end the paper with the main conclu- 
sions and some discussions. 

2. Spherical Collapse Model in QCDM Cos- 
mologies 

2.1. The Basic Equation and Its Solution 

There is a simple way (Weinberg 1972) to derive 
a Friedmann-like equation to describe the evolution 
of the radius of the clusters. Consider a test particle 
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Fig. 1. — Consider a test particle on the edge of a clus- 
ter. Left panel, if dark energy cluster synchronously 
with ordinary matters, then energy conservation law 
will give us a standard Friedman equation eq(l). Right 
panel, if on the scale of galaxy clusters, dark energy 
does not cluster at all, then energy conservation law 
will also give us a Friedmann-like equation, with the 
potential energy of the test particle changed as in 
eq(2). 



It is worth noting that, besides the physical radius 
of the clusters, the symbol a p appearing in eq(l) can 
also be understood as the scale factor of the clusters- 
inside space-time. But that in eq(2) it can only be 
understood as the radius of the clusters, it cannot 
be understood as the scale factor of the cluster-inside 
space-time. This is because, when dark energy is not 
assumed to move synchronously with ordinary matters 
on the galaxy clusters scale, the metric function U (t, r) 
and V(t, r) in the ansaltz (Weinberg 1972) ds 2 = 
-dt 2 + U(t, r)dr 2 + V(t, r)(d0 2 + sm 2 6d(f> 2 ) cannot be 
factorized U(t,r) = a 2 (t)f(r) and V(f,r) = a 2 (t)r 2 
uniformly inside the clusters (Zeng & Gao 2005). 

Combining eq(2) with the background universe 
Friedmann equation: 
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we get 



which moves synchronously with the the edge of the 
clusters at the gravitation of clusters, see FIG.l. De- 
note the scale factor of the background cosmology as o, 
while the physical radius of the clusters as a p . At very 
early times, we can set a p w a. If the physical radios is 
taken to be large enough so that it is equal to the nub- 
ble length, then a p = a throughout the total evolution 
process. If dark energy clusters synchronously with 
the ordinary matters, just as we assumed in (Zeng & 
Gao 2005) then according to energy conservation we 
have 
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which leads to the usual Friedman equation directly. 
If dark energy does not cluster on the scale of galaxy 
clusters but only do so on Hubble scales, just as (Wang 
& Stein hardt 1998) and (Weinber & Kamionkowski 
2002) assumed, then 
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which can also be written as a Friedman-like form, 
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come eq(l) in the limit that the radius of the clusters 
becomes the hubble radius. 



where, as a result of [jgj- \ x =\ = 0, 
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Eq(6) is a first order non-linear ordinary differential 
equation which contains a characteristic parameter £ 
but satisfies a two-boundary condition: 
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This kind of equation can be easily solved by the nu- 
merical method described in (Press et al 1992), chap- 
ter 17. We display our numerical results in FIG. 2. 
Which is also depicted in the figure is the results of 
C(w, Qmb.ta) when dark energy is assumed to move 
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Fig. 2. — £'s dependence on to and f2m6,ta- Left: dark 
energy is assumed only to move synchronously with 
ordinary matters on Hubble scales, on the galaxy clus- 
ter scales, not to cluster at all. Right: dark energy is 
assumed to move synchronously on both Hubble scales 
and galaxy cluster scales. 



synchronously with ordinary matters both on Hubble 
scale and galaxy cluster scales. The actual fact should 
lie between this two extreme cases. But the result of 
(Wang & Stein hardt 1998) does not lie between this 
two extreme cases, so it cannot be thought as an ap- 
proximation of the actual facts. 

Eq(6) can also be solved analytically. Note 



y 

X 



Op O-ta 
d a p .ta 



-C 3 



so 



Let 



y = (ix [l - a(a ta x) + f3{a ta x) 2 + 



(9) 



(10) 



(11) 



substituting eq(ll) into (6) and equating the two sides 
of the resulting equation order by order in x, we can 
get 
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and a similar expression for (3, ... etc. In this paper we 
only need to know a. Some people may argue that if 
progressional solution of eq(6) is not of the form as we 
wrote in eq(ll) or does not exist at all, then our ansaltz 
eq(ll) will be invalid. We would like to point out that 
if such things occur, then when eq(ll) is substituted 
into eq(6) we can not get a self-consistent equation 
with the two sides equated order by order in x. 



2.2. Extracting The Key Parameters of Press- 
Scheter Theory from Spherical Collapse 
Model 

In the ideal model, if there is an over-dense region 
in a flat background universe, then at very early times, 
this region will expand as the background universe 
expands; but because this region's over-dense, its ex- 
panding rate will decrease and stop doing so at some 
middle times; then it starts to shrink because of self- 
gravitating, the final fate of this over-dense region is 
a singular point. But in practice, when this region 
shrinks to some degree, the pressures originate from 
the random moving of particles inside the over-dense 
region will balance the self-gravitation and the system 
will enter the virialization period. In theoretical stud- 
ies, it is usually assumed that the virialization point is 
coincident with the collapse point of the ideal model 
on the time axis. 

According to Press-Shceter theory, if an over-dense 
region is to be virializated at some time a c , its density- 
contrast should be no less than 5 c (w, f2 m Qj a c)- 
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where p mc and p m b are the matter densities of clus- 
ter and background respectively, while Di(a) is the 
growth function of linear perturbation theory (Dodel- 
son 2003), 
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To reduce the numerical computation burdens, we 
will use the fitting formulaes provided by (Ma et al 
1999). We check that when -1.7 < w < -0.4 and 
0.1 < n m o < 0.7, the formulaes provided by this work 
is accurate to 2%. To relate the mass of a galaxy clus- 
ter with its characteristic X-ray temperature, the ratio 
of cluster/background matter densities at the virializa- 
tion point is another very important parameter, 
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It can be shown that Di(a) a ^Q — > a. Using eq(ll), 
we have 
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Substituting eq(16) into eq(13) and using the fact that 
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About A^,s calculation, using definition eq(15) we 
have 
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To calculate the second factor of the above equations' 
right-most part, we can use energy conserving condi- 
tion and virial theorem. If Quintessence clusters syn- 
chronously with ordinary matters, by assuming that at 
the collapse point, the system virializes fully, we can 
get the following relations: 

Ekinetic,c — ~^Umm,c + U m Q :C T^QQ^c (19) 



+ 2J7 m Q jC + - Uqq. c 

= U mmyt a + U m Q^ ta + UqQ.ta (20) 



just as we do in (Zeng & Gao 2005). Now, since 
Quintessence is assumed not to cluster on the scale 
of galaxy clusters at all, we can only write down the 
following relations: 
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in the above four equations, U mm , U m Q and Uqq 
denote the matter-matter, matter-Quintessence and 
Quintessence-Quintessence gravitation potentials re- 
spectively. The subscripts jC and t ta indicate that 
quantities carrying them should take values at the col- 



lapse and turn around time respectively. Explicitly, 
3 GM 2 3 



U„ 



u, 



5 a p . c K 

3 GM 2 3 , * 



ram.ta — 



r _ ^ r^rnc.ta u p,ta 

« u p,ta ° 



c \5 



Uqq.c oc --p Qb ^{a pM — ) 



a ta 



(23) 



UQQ.ta OC -TPQb,ta a l,ta 



U m Q,c OC -[-( 



3 , Clp.ta a c \ 2 5 
1 ~ Tn\Pmc,cPQb, c a p:C 



2 Qp.c ®ta 



10 J 



tt [3,ftp,iTOi a ta \2 ^ l 5 

UmQ,ta OC — [-(^ J - —\Pmc,taPQb,ta a p,ta 

^ ®p,ta Q*ini 



By our normalization, a Pi , n , w Substituting 
these expressions into eq(22) and dividing the result- 
ing equation by the relation p 2 mc ^tf } c = p 2 mc ,t a a% M we 
get: 
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Using energy conservation law and the approximate 
mass conserving condition, p mc ^ c = /w :ta [ v^] 3 , 

Pmc,ta = (PmbM and p Qb , c = PQb,ta [ff ] ~ 3< > 1+w ) , WC 

can change eq(24) into the following form: 
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where £ = |[^] 3 ( 1+M '). Eq(25) can be solved analyt- 
ically, 
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Fig. 3. — For too small f2 m o, if a cluster is to virialize 
too later and w < — 1, then at the virialization point, 
its radius will be larger than that at the turn around 
time, i.e. < 1. 

Using the fact that t c — 2t ta and background Fried- 
man equation (^) 2 oc (p m b + PQb) we can set up an 
integration equation 

r "da' J- ^V^=2 [""da'*/- 

(27) 

where vq = ■ Solve eq(27) numerically, we can 

get the relation at a v.s. a c . Substituting eq(26) and 
— solved from eq(27) into eq(18), we will finally get 
the quantity A c . 

We note that, if w < — 1, then for very small fi m o> 
if a cluster is to virialize too later, then at the virial- 
ization point, its radius will be larger than that of the 
turn around time, i.e. < 1, see FIG. 3. Physi- 

cally this means that, after the "turn around" point, 
instead of collapsing, the cluster experience a period 
of expansion to reach virialization status. Mathemat- 
ically this only means that the turn around point is 
a local minimum instead of a local maximum of the 
cluster radius and has no problem in principle. This 
may be new structure which has not been discovered 
previously, we call it "phantom hole" and leave the 
detailed discussion of this structure for future works. 
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Fig. 4. — S c 's dependence on w, fl m Q and a c . If a c — > 0, 
5 C — > 1.686 asymptotically whatever w and Q m o is. 
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Fig. 5. — A c 's dependence on w, £l m o and a c . As 
a c — > 0, A c — > 178 asymptotically. 

The formation of " phantom hole" will make the kinetic 
energy of the matter-system in it less than 0, please 
see eq(30), which may be a serious problem. 

We provide numerical results for S c (w, fl m o, o. c ) and 
A c (w, fi m o, a c ) in FIG. 4-5. About this two figure, 
what we would like point out is that, as a c — * 0, 
S c — > 1.686 and A c — * 178 whatever w and O m o is. 
Physically, this is because, the earlier an over-dense 
region collapse, the more the background universe is 
like a totally matter dominated one. While in a totally 
matter dominated universe, the fact that 5 C — 1.686 
and A c = 178 can be proved analytically. 



3. The Number Density of Galaxy Clusters 
and Its Evolutions 

3.1. Theoretical Formulaes 

According to Press-Schechter theory, the comoving 
number density of clusters which have collapsed (i.e., 
virialized) at certain red-shift z and have masses in the 
range M <~ M + dM could be calculated (Dodclson 
2003): 
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where: p m b is the matter density of background uni- 
verse; the factor ^jj- denotes the average number den- 
sity of clusters with mass M — ^-R 3 p m b; S c is given 
by eq(13); and 
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In this paper we use the notations of (Dodelson 2003) 
where n s is the primordial power spectrum index, T{k) 
is the BBKS transfer function, W(k ■ R) is the top-hat 
window function. 

To relate the mass of a cluster with its character- 
istic X-ray temperature, consider a virialized spheri- 
cal over-dense region in the background universe con- 
taining dark energies. If dark energies cluster syn- 
chronously with ordinary matters, just as we assumed 
in (Zeng & Gao 2005), then we have 
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where V 2 ir is the mean square velocity of particles in 
the cluster when the system is fully virialized and a p 
is the scale factor of the cluster. However, if dark 
energies do not cluster on the scale of galaxy clusters 
at all, then what we can get should be 
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where a p should be understood as the radius of clusters 
in units of Hubble length if -1 , since in this case we 
cannot define a scale factor globally in the clusters and 
we normalize when a — > 0. Using relation: 
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where ks is the Boltzmann constant, m p is the mass 
of proton, while \xm p is the average mass of particles, 
(3 is the ratio of kinetic energy to temperature. The 
composition ^ has physical meaning of energy trans- 
formation efficiency from thermal dynamic form to x- 
ray form. Substituting eq(30) into eq(31) we get the 
following mass-temperature relation: 
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and A c given by eq(18) and z — a c 1 — 1. 

Just as (Wang & Stein hardt 1998) pointed out, 
since the mass-temperature relation is red-shift de- 
pendent, simply substituting eq(32) into eq(28) cannot 
give us correct number density of clusters in a given 
temperature range today. Instead, we should first find 
out the virialization rate and multiply it by the mass- 
temperature relation then integrate over red-shift 
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From eqs(34), (32) and (29) we can see that be- 
sides j|, n(T,z) will also depend on the cosmological 
parameters w, Cl m o, h, n s and the normalization ag of 
the cosmic density fluctuations. In principle, if we can 
measure the number density v.s. temperature relation 
precisely enough, by numerical fittings, we can deter- 
mine all these parameters simultaneously from obser- 
vations. However, in practice, because of parameter 
degeneracy and measure errors, we can only determine 
some of them or their combinations partly. 

Now let us return to the strange phenomenon dis- 
played in FIG.2 and 3, i.e. C < 1 and < 1 re- 
spectively. We have explained that these two things 
take place when w < — 1 and £l m b,ta or Q m o takes too 
small values. It can be checked that when < 1, 
we must have ( < 1 and the kinetic energy of the viri- 
alizcd matter system will be less than zero, see eq(30). 
We think such a "cluster" can not emit X-rays. 

3.2. Numerical Results, Effects of w on the 
Number Density of Galaxy Clusters and 
Its Evolutions 

We display the effects of w on the number-density 
v.s. temperature of galaxy clusters in FIG. 6 for four 
compositions of cr 8 and ^. For any given parameter 
set {cr 8 , p}, if — 1 < w < 0, then w affects the number- 
density v.s. temperature of galaxy clusters exponen- 
tially. If w < — 1, the effects are weak. 

In FIG. 7 we display the effects of w on the number- 
density v.s. red-shift of galaxy clusters whose mass is 
greater than 8/i _1 x 10 14 M sun for two values of cr 8 . 
From the figure we can easily see that w affects the 
number-density v.s. red-shift relation of galaxy clus- 
ters remarkably. When — 1 < w < 0, the number 
density almost does not vary with time for <r 8 = 0.55, 
but for cr 8 = 0.85 case, the number density increases 
as we look back to the past. When w is less than 
— 1, the number density decreases at low red-shift, but 
increases at high red-shift, the turn around red-shift 
depends on w. The less is w, at the higher red-shift 
the trend reverses. In FIG. 8, we depict the same ef- 
fects for the number density of galaxy clusters whose 
mass is greater than 1.5/i _1 x 10 14 M stm . From the fig- 
ure, we can see for the smaller mass galaxy clusters, 
the effects of w on the number-density v.s. red-shift is 
more remarkable. 

Comparing this fact with the observational results 
reported by (Bahcall, Fan & Cen 1997) and (Bahcall 
& Bode 2003), we can almost immediately say that w 



can not be greater than — l! Of course, a reliable con- 
clusion should be obtained by best fitting the current 
observational results with theoretical predictions. But 
if we want to fit the results of (Edge et al 1990; Hen- 
gry & Arnaud 1991; Henry 1997) with our theoretical 
formulae eq(34), we at least have four parameters ^, 
w, cr 8 and fl m o to determine. If we want to fit the 
results of (Bahcall & Bode 2003), we have to treat 
the problem of fitting data with errors in both coor- 
dinates. This two kinds of operations both cost time 
formidably. We leave them for future works. 

Comparing the results in this paper with that of 
(Zcng & Gao 2005), we can see that the effects of w 
on the evolution of number density of galaxy clusters 
is even more strong under the assumption that dark 
energy does not cluster on the scale of galaxy clusters 
than the case where dark energy is assumed to clus- 
ter synchronously with ordinary matters. It's easy to 
imagine that, the actual case should lie between this 
two extreme way, matter clusters and forms potential 
well, then dark energy falls into it and cannot climb 
up so clusters also. If we want to fit observation re- 
sults into theoretical predictions to get reasonable con- 
straints on w, f2 m o, cr 8 and other cosmological param- 
eters, theoretical formulaes in the both extreme cases 
are needed, and some kind of interpolation should be 
used. 

Whatever the actual case is, according to the results 
found in this paper and its previous counterpart (Zeng 
& Gao 2005), we know that the effect of w on the 
number density of galaxy clusters is so strong that we 
think it should be possible to use this effect to measure 
w. 

4. Conclusions 

We study the top-hat spherical collapse model 
of galaxy clusters formation in the flat QCDM or 
Phantom-CDM cosmologies under the assumption 
that Quintessence or Phantom docs not cluster on 
this scale. We find that under this assumption, the 
key parameters of the model exhibit rather non-trivial 
and remarkable dependence on the equation of state 
coefficients w of Quintessence or Phantoms. We then 
applied the results in Press-Scheter theory and calcu- 
lated the number density v.s. temperature function 
and the evolution of the number density of massive 
galaxy clusters and find that these two Quantities are 
both affected by w exponentially. 

For the number density v.s. temperature function of 
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Fig. 6. — Effects of w on the number density v.s. tem- 
perature function of galaxy clusters when z = 0. The 
larger is a$ or 4, the larger the function value will 
be. All four figures have f2 m o = 0.27, h — 0.71 and 
n s = 1.0. 

galaxy clusters, we find that it is an increasing func- 
tion of w and the dependence on w is more strong 
in the range — 1 < w < than it is in the range 
— 00 < w < — 1. While for the evolution of the num- 
ber density of massive galaxy clusters, we find that 
when w is less than —1, the number density decreases 
at low red-shift, but increases at high red-shift, the 
turn around red-shift depends on w. The smaller is 
the galaxy clusters' mass, the stronger is the this ef- 
fect. According to the observational result, which says 
that the number density of massive galaxy clusters de- 
creases as we look back to the past, we can qualita- 
tively conclude that w should not take values greater 
than —1 too much. It should take values less than — 1. 

The actual dark energy cluster property should lie 
between the two extreme cases discussed in this pa- 
per and that in (Zeng & Gao 2005). But whatever 
the fact is, our results here and that in (Zeng & Gao 
2005) indicate that, w affects the number density of 
galaxy clusters exponentially So we think measuring 
the number density of galaxy clusters and its evolu- 
tions may be an effective way to determine w. 

As discussions, we would like to state that, if the 
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Fig. 7. — Effects of w on the evolution of number 
density of galaxy clusters whose mass is greater than 
8/i _1 x 10 14 M sun . When w is greater than — 1, the 
number density almost does not vary with time for 
erg = 0.55, but for erg = 0.85 case, the number den- 
sity increases as we look back to the past. When w 
is less than —1, the number density decreases at low 
red-shift, but increases at high red-shift. Both figure 
have fi m0 = 0.27, h — 0.71 and n s = 1.0. 

l°g 10 n M>1.5 tr'xlO^Wi lo Sl0 n M>1.5 tr'xlo'^M™ 



o-s=0.55 o-8=0.85 




Fig. 8. — The same as FIG. 7, but we consider the num- 
ber density of galaxy clusters whose mass is greater 
than 1.5/iT 1 x 10 14 M SU „. 

problem we pointed out in (Zeng & Gao 2005) is the 
fact, then we may have to accept that, we ignored a 
very important assumption in the usual ACDM cos- 
mology. That is: the cosmic component denoted by A 
and leading to the accleration of the universe moves 
synchronously with ordinary matters. The reason is 
very clear, if A does not move synchronously with or- 
dinary matters on Hubble scales, then in our comoving 
reference frame build on the ordinary matters(such as 
supernovaes) , we should have a A current flowing out- 
side the hubble horizon. Once that occurs, we cannot 
define an all universe uniformly defined scale factor, so 
we have no Friedmann equation at all. In that case, 
our current explanation of the acceleration indicated 
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by the observation of supernovaes will be very prob- 
lematic. 

If A component moves synchronously with ordinary 
matters, then we will have no reason to say that it is 
the vacuum energy of quantum field. On the contrary, 
some kinds of couplings between A (or dark energy) 
and the ordinary matters is a must be derivation. 

Part of the numerical computations are performed 
on the parallel computers of the Inter-discipline Center 
of Theoretical Studies of ITP, CAS, Beijing, China. 

Appendix. Comparing With Previous Works 

A. Comparing The Basic Equations 

Our results in this paper and its sibling one (Zcng 
& Gao 2005) are so different from the previous works 
that without a concrete comparison and explicit point- 
ing out the problem in those works, few peoples will 
believe our conclusions. 

To compare our basic equations with the basic equa- 
tions used by (Wang & Stein hardt 1998), we can dif- 
ferentiate eq(2) and use energy conservation law to 
get, 



^ = -4ttG 



1 3 a a 1 {pQ b a 2 p y 1 

(w + dPQ&~- + a • + Q^ mc 

A 2 a p a p 6 a p a p 3 



In the right hand side of this equation, if we let 
a == a pi i.e., assume that the dark energy moves syn- 
chronously with ordinary matters on the galaxy cluster 
scales, then it reduce to the eq(A2) of (Wang & Stein 
hardt 1998). Note, in the right hand side of eq(35), a 
only appears in terms where dark energy is involved. 
But (Wang & Stein hardt 1998) does not use this as- 
sumption consistently, because when it combine its 
eq(A2) with the equation describing the dark energy's 
evolution eq(A6), it uses the assumption that only on 
Hubble scales, dark energy moves synchronously with 
ordinary matters. 

Although we have pointed out in (Zcng & Gao 
2005), we still would like to point out that, as long as 
variable separation technique is used in solving Ein- 
stein equation, then whichever (Raychaudhuri or Frei- 
dmann) equation we choose to describe the evolution 
of the over-dense region, we will in fact have assumed 
that dark energy moves synchronously with ordinary 
matters in our over-dense region. Let us explain this 



point in more details. For an over-dense collapsing re- 
gion, the most general metric describe its inside space- 
time is 



ds 2 



(35) 



-dt 2 + U(t, r)dr 2 + V(t, r)(d9 2 + 



sm 



(36) 



Using Einstein equation G M „ = SirGT^, we can 
prove that, only when no energy current flowing out- 
side the over-dense region, is the energy momentum 
tensor T^ u diagonal, and can the metric function 
U(t,r) and V(t,r) be factorized as 



U(t,r)=a 2 p (t)f(r), V(t,r) = a 2 p (t)r 



(37) 



And only when the U (t, r) and V (t, r) function is fac- 
torized, can we have 



4irG 



[pmc + PQc + 3pQ c ], 



a p 3 
in (WS98)'s notation 

R An/ , PQ + Pcluster s 
- = -^G{ P Q + -= ) 



(38) 



dp a v . „ k 

2 ~ 2 

Qp Qp O'p 



2— =4nG[p mc + pq c -pq c ]. (39) 



As a must be of eqs(38) and (39) 



8ttG 



{Pmc + PQc) 



(40) 



So Frcidmann equation and Raychaudhuri equation 
must hold at the same time, or must not hold simulta- 
neously. The statement of (Wang & Stein hardt 1998) 
and (Weinbcr & Kamionkowski 2002) that when dark 
energy does not cluster on the scale of galaxy clus- 
ters, Raychaudhuri equation can be used to describe 
the evolution of the over-dense region but Freidmann 
equation does not hold is an incorrect statement. 

Let us say more explicitly, when dark energy does 
not cluster on the galaxy clusters, hence a dark energy 
exists which describe the flowing of dark energy out- 
side the over-dense region, the basic equation which 
should be used to describe the evolution of the over- 
dense region is not the cq(A2) of (Wang & Stein hardt 
1998), it should be our eq(2). Our eq(2) is not ob- 
tained by variable separation in solving Einstein equa- 
tion. We obtained it by energy conservation, so we 
include the effects of the dark energy current on the 
evolution of the over-dense region. 
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B. Comparing The Numerical Results 

Let us emphasize again that the problem in the ex- 
isting works (Wang & Stein hardt 1998) and (Weinber 
& Kamionkowski 2002) is, when writing down the ba- 
sic equations describing the evolution of the over-dense 
region they assumed that dark energy moves syn- 
chronously with ordinary matters, please see eq(A2) 
of (Wang & Stein hardt 1998), but when writing down 
equations which will be used to describe the density 
of dark energies in the over-dense regions, please see 
eq(A6) of (Wang & Stein hardt 1998), a different as- 
sumption is made. That is, dark energy only moves 
synchronously with ordinary matters on Hubble scales. 

Just as we pointed out in (Zeng & Gao 2005) and in 
the conclusion section of this paper. In realities, dark 
energy should have some degree of cluster behaviors 
on the galaxy clusters. We can imagine, ordinary mat- 
ters cluster and form potential wells, when dark energy 
falls in and some degree of dark energy's clustering will 
occur either. So the actual case of dark energy's clus- 
tering phenomenon should lie between the following 
two extreme cases. The two extreme cases are, dark 
energy moves synchronously with ordinary matters on 
both Hubble scales and galaxy cluster scales or dark 
energy only moves synchronously with ordinary mat- 
ters on Hubble scale but could not fall in the potential 
wells formed by the over-dense matter region at all. 

We study the first extreme case in (Zeng & Gao 
2005) and the second extreme case in this paper. A 
natural question is, can the results of (Wang & Stein 
hardt 1998) lie between our two extreme cases? If this 
is the case, then although inconsistent, the treatment 
of (Wang & Stein hardt 1998) can be thought as some 
kinds of approximation of realities. We will see in the 
following that this is not the case. We compared the 
results of ('s dependence on w and $l m b,ta in FIG. 9 and 
10, FIG. 9 is 3-dimensional, FIG. 10 is 2-dimensional. 

C. Some Comments On References 

When this paper and its sibling one (Zeng & Gao 
2005) are put on the e-preprint arXive and submit- 
ted to Astrophysical Journals, we are told that just 
recently, many authors have studied this problem in 
different depth. Such as (van de Bruck & Mota 2005), 
(Battye & Weller 2003), (Koivisto 2005), (Mota & van 
de Bruck 2004), (Nunes & Mota 2004) and (Manera 
& Mota 2005). We believe there must be more we do 
not know at this moment. We think we should first 



<T(w,a 




Fig. 9. — C' s dependence on w and £l m b,ta- The upper, 
solved from eq(A2) and (A6) of (Wang & Stein hardt 
1998), when the basic equation (A2) were written to 
describe the evolution of the over-dense regions, the as- 
sumption that dark energy moves synchronously with 
ordinary matters was made, but when the dark energy 
density was involved in the right hand side of the equa- 
tion, a different assumption is made. That is, only 
on Hubble scales, dark energy moves synchronously 
with ordinary matters. Down left, dark energy is as- 
sumed to move synchronously with ordinary matters 
both on Hubble scales and on galaxy clusters scale. 
Down right, dark energy is only assumed moves syn- 
chronously with ordinary matters on Hubble scales, it 
is not assumed to move synchronously with ordinary 
matters on galaxy clusters scale. 



express our thanks to the referee and editors of APJ 
and the authors of these papers for their informing us 
of these works, we then also would like to point out 
that, almost all these works contain the problem we 
pointed out in this paper, i.e., when writing down the 
basic equations to describe the evolution of the over- 
dense regions, the assumption that dark energy moves 
synchronously with ordinary matters is made, but in 
the basic equations when the dark energy's density is 
involved, the results from another different assumption 
is used. 
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Fig. 10. — The same as FIG. 9, but fl m b,ta is set to 
three special value 0.1, 0.3, 0.5. 
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